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Abstract 

We investigate existence and properties of trapped surfaces in two 
models of collapsing null dust shells within the Gibbons-Penrose construc- 
tion. In the first model, the shell is initially a prolate spheroid, and the 
resulting singularity forms at the ends first (relative to a natural time 
slicing by flat hyperplanes) , in analogy with behavior found in certain 
prolate collapse examples considered by Shapiro and Teukolsky. We give 
an explicit example in which trapped surfaces are present on the shell, 
but none exist prior to the last flat slice, thereby explicitly showing that 
the absence of trapped surfaces on a particular, natural slicing does not 
imply an absence of trapped surfaces in the spacetime. We then examine 
a model considered by Barrabes, Israel and Letelier (BIL) of a cylindrical 
shell of mass M and length L, with hemispherical endcaps of mass m. 
We obtain a "phase diagram" for the presence of trapped surfaces on the 
shell with respect to essential parameters A = M/ L and ^ = m/M. It is 
found that no trapped surfaces are present on the shell when X or jj, are 
sufficiently small. (We are able only to search for trapped surfaces lying 
on the shell itself.) In the limit A ^ 0, the existence or nonexistence of 
trapped surfaces lying within the shell is seen to be in remarkably good 
accord with the hoop conjecture. 



1 Introduction 

It is well known from the singularity theorems Q that singularities must occur in 
a wide variety of circumstances relevant to gravitational collapse — in particular, 
when trapped surfaces or reconverging light cones are present. However, there 
is still little direct evidence on whether cosmic censorship holds, i.e., whether the 
singularities of gravitational collapse are always hidden in black holes (see |^ and 
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[|| for recent reviews). Much of our knowledge of what occurs in gravitational 
collapse is based upon the analysis of spherically symmetric models and their 
linearized perturbations. Without results from other models, we cannot know 
which features might be shared by — and which absent from — general kinds of 
collapse. 

It has been suggested by a number of authors that highly prolate gravi- 
tational collapse may be of a qualitatively different character from spherical 
collapse (see, e.g. [||). An explicit statement of this view has been expressed in 
the "hoop conjecture" ||^: 

Horizons form when and only when a mass M gets compacted into 
a region whose circumference in EVERY direction is C ^ AttM. 

The intended interpretation of this statement is that a distribution of matter 
that is elongated in one dimension but is sufficiently narrow in the other dimen- 
sions might be expected to collapse without the creation of a event horizon, so 
that any singularity that formed would be "naked" . However, there are serious 
difficulties in giving a precise statement of the hoop conjecture, both in defining 
the "mass" — one would want to include concentrations of gravitational energy, 
for which there is no local definition — and the "circumference" — which always 
can be made arbitrarily small if no choice of spacelike slicing is specified. 

Nevertheless, numerical simulations by Shapiro and Teukolsky ^ of the 
collapse of a collisionless gas cloud appeared to lend support to the hoop con- 
jecture. For a sufficiently prolate cloud, a singularity was found to form (after 
which point the numerical evolution could not be continued) but no trapped sur- 
faces were found prior to the singularity on the family of maximal hypersurfaces 
they used for their time evolution. 

Barrabes, Israel and Letelier considered a simple analytic counterpart of 
the Shapiro- Teukolsky examples within the Gibbons-Penrose construction (see 
below). This model consists of an imploding shell of photons (null dust) in the 
shape of a cylinder with hemispherical cndcaps and with pieccwisc uniform mass- 
energy density. They found, similarly, that given a fixed amount of matter, if 
the cylindrical portion of the shell is made long enough (i.e. sufficiently prolate), 
the 2-surfaces formed by the intersection of the shell with flat hyperplanes of 
the Minkowski region interior to the shell are never trapped at any stage of 
collapse. 

However, in both the numerical and analytic models, the conclusion that no 
trapped surfaces exist relies on a specific choice of time slicing. Indeed, as a 
counterargument to the Shapiro- Teukolsky result, Wald and Iyer |^ proved that 
even in Schwarzschild spacetime, it is possible to choose a (highly non-spherical) 
time slice which comes arbitrarily close to the singularity, yet for which no 
trapped surfaces are found to its past. This example explicitly demonstrates 
that the absence of trapped surfaces in a particular time slicing does not estab- 
lish the absence of trapped surfaces in the spacetime. Nevertheless, the slicing 
in is rather contrived, and it is far from clear that anything similar to this is 
occurring in the Shapiro- Teukolsky examples. 
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The main purpose of this paper is to investigate the existence and proper- 
ties of trapped surfaces in highly prolate collapse in the context of spacetime 
models produced by the Gibbons-Penrose construction. The Gibbons-Penrose 
spacetimes contain imploding shells of null dust, with the region of spacetime 
interior to the shell being flat. These models have the advantage that it is 
relatively easy to search for trapped surfaces lying on (or, more precisely, in- 
finitesimally outside of) the collapsing shell. Unfortunately, in these models, the 
spacetime outside of the shell is not explicitly known, so one cannot determine 
the existence or properties of trapped surfaces lying at a finite distance outside 
of the shell. 

In section 2, we review the Gibbons-Penrose construction. In section 3, 
we examine the collapse of an initially prolate spheroidal shell. We obtain 
an explicit example — previously suggested in ^ — in which, in the natural flat 
hyperplane time slicing in the region interior to the shell, the singularity forms 
at the poles before any trapped surface has finished forming along the equator. 
Thus, a numerical code using this natural time coordinate — which would halt 
at the time the singularity begins to form — would not find any trapped surfaces 
on the shell. Nevertheless, in the fully evolved spacetime, trapped surfaces do 
occur on the shell. Thus, this example explicitly displays the behavior closely 
analogous to that constructed in , but does so in the context of prolate collapse 
with a natural choice of time slicing. 

In section 4, we re-examine the Barrabes-Israel-Letelier (BIL) model of a 
collapsing cylindrical shell with hemispherical endcaps, in the case where A = 
M/L < 1/8, where M is the mass of the cylinder and L its length. BIL showed 
that in this case, none of the 2-surfaces formed by the intersection of the shell 
with the fiat hyperplanes of the interior region are trapped surfaces. Here, we 
search for trapped surfaces on the shell which do not lie in these hyperplanes. We 
present a "phase diagram" for the the presence or absence of trapped surfaces 
as a function of the essential scale-invariant parameters A and = m/M, where 
m the mass of each endcap. We find that no trapped surfaces are present on 
the shell for A or /i sufficiently small (although we expect that trapped surfaces 
will be present outside of the shell in all cases). We find that the existence 
or nonexistence of trapped surfaces lying on the shell is in remarkably good 
agreement with what would be expected from the hoop conjecture. 

2 The Gibbons-Penrose construction 

The Gibbons-Penrose construction considers a thin shell of matter collapsing 
inward at the speed of light from past null infinity. The trajectory of the shell 
describes a null hypersurface Af which divides the spacetime into two regions: 
the interior, which is Minkowski, and the exterior, which is not flat and generally 
contains gravitational radiation. Singularities form where the null hypersurface 
develops caustics. 

The construction is completely specified by two functions, one defining the 
shape of the collapsing shell, and one defining its (distributional) mass density. 
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The degrees of freedom associated with the shape correspond to a choice of 2- 
surface S(T) in the hyperplane t — T of Minkowski spacetime. The hypersurface 
J\f is then generated by the ingoing normal nuU geodesies through S{T). We 
restrict to the case where S{T) is convex, so that no caustics appear to the past 
of our initial time. The "initial" mass density (t(T) can be chosen to be any 
non- negative function on S{T). The mass density is then extended to the whole 
of M by the conservation equation. 

One can then, in principle, go about searching for trapped surfaces in the 
spacetime. In general, this, of course, requires solving Einstein's equations in 
the region of interest. The great advantage of the Gibbons-Penrose construction 
is that as long as one focuses attention on surfaces infinitesimally outside A/", 
one need not know anything about the spacetime outside A/"; that is, given some 
candidate compact 2-surface in A/", we can determine whether it is trapped 
without ever solving for the exterior spacetime. The expansion 9 of future- 
directed, outgoing null geodesies normal to the 2-surface, just before they cross 
the shell, is @ 

6'insidc = -Tr L (1) 

where L is the extrinsic curvature of the 2-surfacc under its embedding in 
Minkowksi spacetime, with respect to outgoing null normals. Just outside the 
shell, we have 

"outside — t^insidc matter (2) 

where, by integrating the Raychaudhuri equation jl^ across the shell (so that 
only the distributional source term contributes), we find that 

(A0)„,attor = -87ra (3) 
Thus, the 2-surface is trapped if 

- Tr L < Sttct (4) 

everywhere on the surface. 

The project is further simplified if one looks instead for an outer marginally- 
trapped surface (MTS), i.e. a compact 2-surface for which the expansion of 
future-directed outgoing null normals is everywhere zero. By an argument sim- 
ilar to one given in [Q for trapped surfaces lying in a spacelike hypersurface, if 
M contains trapped surfaces, then the outer boundary of the region containing 
the trapped surfaces is an MTS. The absence of an outer marginally-trapped 
surface on M thus implies the absence of trapped surfaces on M . 

The advantage of searching for MTS's is that a MTS satisfies the equation 

- Tr L = Sttct (5) 

rather than an inequality. Thus, rather than attempt to explore the myriad 
of possibilities for where trapped surfaces might occur, one studies eq. and 
attempts to determine if any global solutions exist. 

In the following, we narrow our search for trapped surfaces in the whole 
spacetime to the search for MTS's confined to the null hypersurface N . 
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3 The collapsing spheroid 



We now consider the collapse of a prolate ellipsoid in the Gibbons-Penrose 
construction. We choose Minkowski coordinates {t, x, y, z) in what will be the 
interior of M , and specify the initial shape iSq at t = to be an axisymmetric 
ellipsoid (i.e., a spheroid) 

x'^{9,(j)) — (0, a sin 6 cos cf), a sin 6 sin 0, b cos 9) (6) 

e [O,7r],0e [0,27r) (7) 

where a < 6 so that the spheroid is prolate. Following the ingoing null normals 
forward and backward in time yields the form of A/", 

x'^{t, 6, (j)) — (t, (a — "^^"j '^i'^ ^ ^'^'^ 4*' 0^ ~ sin 6* sin 0. ^& — '~r^ '^'^^ ' 



(8) 

where 

E{0) ^ cos^ 9 + b^ sin^ e (9) 
The metric induced on the intersection of Af with the hyperplanes of constant 

t is 

abt \ „ ,„o 2 f bt \ .2 



-[^~^) Ede^ + « (^1 - ^ j (10) 

This metric becomes truly singular when the quantities in parentheses become 
zero. This is just the condition that t is equal to a principal radius of curvature 
of the initial spheroid, where these radii of curvature are given by 

^r^--i(^)-^ (11) 

= = JtJ2 (12) 

Since Ki < K2 (with equality only at the poles), the singularity actually occurs 
at t{e) = l/K2{e). 

The evolution of the initial mass density ctq (which so far is unspecified) is 
governed by the local mass-conservation equation 

|(adA)=0 (13) 



where dA is the area form obtained from the metric (10), 

dA = aUVVE sin 9 d9 A dcj) (14) 

where 

U{9,t) = l~K2t (15) 
V{9,t) = l-Kit (16) 
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Given some ao, then, 

a{t,e,<p)^^ (17) 

Equations and define the Gibbons-Penrose construction for the collaps- 
ing spheroid. 

Now consider an arbitrary spacelike 2-surface on J\f, 

t^F{9,(t>) (18) 

We seek the condition on F such that this surface is marginally-trapped. The 
embedding of the 2-surface in Minkowski spacetime is given by equations ( p^ 
and ( p|), and its extrinsic curvature is readily calculated. Combining this result 
with (|17|), we find that eq. becomes 



where Laplacians and gradients are with respect to the metric and the 
functions U and V are to be evaluated on the 2-surface {i.e. replace t by F{9, (p) 
everywhere) . 

Finally, we will restrict attention to axisymmetric solutions F = F{9) since 
the outer boundary of the region containing trapped surfaces must have the 
same symmetries as the spacetime itself. In this case eq. (|l^) reduces to 

U d /sme\ ,„ „ _ FF' f F' E' 



\ FF' fF' E'\ 



+ ^(Ki-l-K2-2KiK2F) = 87rcro (20) 



The parameters and free functions of this model can be made to yield a 
spacetime where the outermost marginally-trapped surface encloses the singu- 
larity, but forms later at the equator than the singularity does at the poles, in 
the time coordinate t. A numerical code using this time coordinate would then 
halt before any trapped surfaces have formed. To obtain an explicit example 
where this occurs, we will specify a form of i^, and use ( |20| ) to determine the 
requisite initial mass density, which we must verify to be everywhere positive 
and finite. A relatively simple choice of MTS is the surface of constant V 

Fm = ^ (21) 

where 7 is some constant. In order that the MTS have the special properties 
described above, we must restrict 7 according to 

^ < 7 < 1 (22) 

Substituting F into ( pO| ) gives the initial mass density, 

CTo = ^(A cos^ e + B cos^ + C cos^ e + D) (23) 
Stt 
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Figure 1: The chosen MTS (7 = 1/2, a — 1, b — a/5) and the singularity (0- 
coordinate suppressed). Dotted lines indicate null generators of the hypersurface 
Af. Note that the MTS forms later at the equator than the singularity forms at 
the poles. 

where A, B, C and D are constants depending on 7, a and b. For the choice of 
parameters 7 = 1/2, a = 1, b = for example, we find A = —352/225, B = 
8/9, C = 2/9, and D = 11/6, from which it can be verified that uo is everywhere 
positive. The MTS for this case is plotted in fig. 0. 

The possibility exists, however, that our chosen MTS may not be the outer- 
most one, and so would not be the outer boundary of a trapped region. (Exam- 
ples where there exists an MTS which is not the outer boundary of a trapped 
region will be given in the next section.) In that case trapped surfaces would 
occur to the past of our MTS, and our argument about the inability to find 
trapped surfaces would not apply. We have investigated this issue numerically 
and have found that for at least some choices of parameters a, b, and 7, our 
MTS is indeed the outermost one. 

We numerically examined the boundary value problem consisting of ( pO| ) 
together with F'{-k/2) = (to ensure reflection symmetry about the 6* = 7r/2 
plane) and -F'(O) = (to ensure regularity at the poles). The singular nature of 
the differential equation aX 6 — Q prevents the direct application of the boundary 
conditions in the numerical routines. Instead, the equation is solved from 9 = 
7r/2 toward 9 — 0, with initial conditions F'{tt/2) — and F{tt/2) = to. We 
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then check to see if the boundary condition F'{Q) = is met. 

For the choice of parameters 7 = 1/2, a = 1, b = a/S, the chosen MTS has 
to = 0.5. For to slightly less than 0.5, the solutions "curve more to the past" 
and — > cx) as 6* ^ 0. As to is further decreased, the solutions diverge 
even more rapidly as 9 ^ 0, and the boundary condition F'(0) = cannot be 
met. We conclude that the chosen MTS is the outermost one, and therefore no 
trapped surfaces exist outside of our chosen MTS.F] 



4 The Barrabes-Israel-Letelier model 

In the Barrabes-Israel-Letelier model (BIL) model of prolate collapse, the ini- 
tial shape is a cylinder of length L with hemispherical endcaps. As this shell 
collapses, it retains this form; this is one advantage of this model. The resulting 
singularity will be a "spindle" singularity with stronger, point singularities at 
the ends. 

Again, we let (t, x, y, z) be Minkowski coordinates for the region interior to 
the shell. Since the singularity forms "simultaneously" , it is convenient to label 
the hyperplane at which the singularity forms as t = 0, and choose t to increase 
towards the past (so that t is positive in the region of interest). We can write 
down A/" immediately in piecewise form: 

x^(t,z,(j)) = {— t,t cos (j3,tsiTi(f>,z) cylinder 

a;'^(t, 0, 0) = (— i sin0 cos t sinfl sini/), t cos 6') right endcap (24) 

x'^{t,9' — (— sin 0' cos (/>,t sin 0' sine/), f cos 0') left endcap 

^e[-^,|];0,0'e[O,|];0e[O,2^) (25) 

The pieces are matched by identifying z — L/2 with 6 = tt/2, and z = — L/2 
with 9' = 7r/2. 

The initial mass density is taken to be uniform over each component of the 
collapsing shell. Because the shape is retained during collapse, the density varies 
with time as 

We will refer to M as "the total mass-energy of the cylindrical portion" , and m 
as "the total mass-energy of each endcap" . Note that, except for one value of 
t, the mass density is discontinuous where the cylinder meets the endcaps. 

Again, we seek MTS's that have the same symmetries as Af, since the out- 
ermost MTS (which is the one of interest to us) must have these symmetries. 
Because of the reflection-symmetry of the MTS about z = 0, we focus on z > 

^For to slightly larger than 0.5, the solutions "curve more toward the future" and F'{8) — > 
—oo as 9 — > 0; these surfaces can be perturbed slightly to construct regular trapped surfaces 
bounded by our chosen MTS. When to is slightly larger than 0.51, solutions run into the 
singularity before reaching 6 = 0. We therefore also conclude that there is no MTS forming 
at a later time than our chosen MTS. 
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and ignore the left endcap. Then, exploiting the axial symmetry, we can define 
a candidate 2-surface by the equations 



[ V(0) endcap ^ ' 

It has been shown ||^ that if / and V are to describe an MTS, they must satisfy 

dz'^ \dz I L 



+ =^ (28) 



The MTS must be smooth at z = 0, regular at 61 = 0, but need only be at 
z = -J since A/" itself is only and the density is discontinuous there. Hence, 
we require 

|(0)^0 (30) 

v'(f ) = .n^) (31) 

d ,7r, df ,L. . . 

;jj(logF)(-) = -^(-) (32) 

^(0) = (33) 

The model is completely specified by the three parameters m, M and L. 
Since, as m, M and L are scaled together, the solution will scale with them, it 
is useful to define A = MjL and = m/M . Then a given pair (A, ji) defines an 
equivalence class of solutions with respect to overall change of mass/distance 
scale. 

In the original presentation of the BIL model 0, it was demonstrated that 
for A < 1/8, "the collapsing shell" (meaning the intersection oi N with a flat 
hyperplane) is never a trapped or marginally-trapped surface, whereas for A > 
1/8, these intersection surfaces eventually become trapped. Our aim here is to 
investigate whether in the case A < 1/8 trapped surfaces exist on the shell that 
do not lie on flat hyperplanes. We therefore restrict attention to A < 1/8. 

To investigate this issue, we proceed as follows. First, we fix M and L and 
solve eqs. ( p8| ) and ( pO| ) on the interval [0, L/2], with to = /(O) treated as an 
arbitrary free parameter. (Note that the parameter m does not enter these 
equations.) A unique solution to these equations exists |^ for all choices of 
^0 > 0. The resulting value of / and its derivative at z = L/2 then generate the 
required initial conditions (|l|) and (|3|) for the system (|9|) , (|l|) , (^ and (||). 
This is an overdetermined system, so, in general, no solution will exist. However, 
in appendix A, we outline a proof (due to David Ross) that there always exists 
a unique choice of m such that a solution to these equations exists. 
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Figure 2: Qualitative behavior of m vs. to. 



For each fixed M and L, this unique solution defines a continuous function 
mlto). The asymptotic behavior of this function as io ~^ and as to —^ oo is 
analyzed in appendix B. As also noted there, it is not difficult to show that no 
solutions exist for m — 0. This fact together with the continuity and asymptotic 
behavior of m{to) implies that there exists an mo > such that m{to) > mo for 
all to- The resulting behavior of the function m(to) is sketched in fig. ^. 

From inspection of fig. ||, we see that at each fixed M and L, there exist three 
qualitatively distinct regimes with regard to solutions of the MTS equation. For 
< m < mo, there do not exist any MTS's lying on M . On the other hand, for 
Too < to < TOi, there exist at least two MTS's on M . Finally, for to sufficiently 
large, there exists a unique MTS lying on N. 

In order to gain more insight into the behavior of the MTS's, we have nu- 
merically studied the solutions to the MTS equations for a wide range of the 
essential parameters A and /x. For any fixed > 0, we find that when A is suf- 
ficiently close to 1/8, there always exists a unique MTS (see fig. |). This MTS 
is the outer boundary of a trapped region. As A is decreased holding /i fixed 
{i.e., as the shell is made more prolate), one encounters a critical value, A = Ai, 
beyond which a second MTS appears inside the first one, almost touching the 
singularity. As A is decreased further, the outer MTS and inner MTS migrate 
towards each other, until, at a second critical value A = Ao, the two surfaces 
merge and then disappear (see fig. ^). For A < Ao, no marginally-trapped sur- 
faces are found confined to the null hypersurface. These critical values Ao and 
Ai describe curves in the (A, /i) plane, dividing parameter space into regions for 
which qualitatively different types of solutions occur (fig. H). 

The disappearance of trapped surfaces on M when A < Ao (/i) is in accord 
with the qualitative behavior one might have expected from the hoop conjecture. 
To investigate this more quantitatively, we interpret the "mass" appearing in 
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Figure 3: Interpreting the plots. We represent the MTS on M by projecting 
it onto a flat hyperplane in the Minkowski region and then suppressing the (j)- 
coordinate. Shown is a typical MTS (/U = 2,A = 0.08). In the central region 
(to the left of the dotted line), the ordinate yields the time to the past of the 
singularity on null generator labelled by z on the cylindrical portion of the 
surface. In the end region (to the right of the dotted line), the radius (measured 
from the intersection of the dotted line with the abscissa) as a fimction of the 
polar angle 9 yields the time to the past of the singularity on null generator 
labelled by 6 on the hemispherical endcap. 
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Figure 6: The same phase diagram as in fig. with the function A — (27r(l 
also plotted (dashed line). 
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the general statement of the hoop conjecture (see section 1) as being the total 
mass-energy, Mtot = M + 2m, of the shell, and we take the "circumference" of 
a 2-surface on M to be the pole-to-pole circumference 2L -\- 41, where 

7r/2 

V{e)d9 (34) 

Since C > 2L for all 2-surfaces on Af, with equality holding to arbitrarily good 
accuracy as one approaches the singularity, if this interpretation of the hoop 
conjecture is correct, one would expect that trapped surfaces will exist on Af if 
and only if 2L < 47r(M -|- 2m), i.e. if and only if 

i < 27r(l + 2/i) (35) 
A 

Thus, the hoop conjecture suggests that the curve defined by equality in eq. 
( ^5| ) should be the "dividing line" between the regimes where trapped surfaces 
do and do not occur on Af. This curve is plotted as a dashed line in fig. ^. 
It can be seen from the figure that this line lies remarkably close to the actual 
dividing hne^ A = Ao(/i). 

Furthermore, even better quantitative agreement with expectations from the 
hoop conjecture is obtained if we consider the ratio C/Mtot for the "last" MTS, 
i.e. if we examine, as a function of /i, the quantity 

^ = ^^ + ^^ = ^ + ^ (36) 
Mtot M + 2m (1 + 2/i) • Ao(Ai) ^ ^ 

where / is given by eq. (|3^ ) above. In numerical solutions, we find that the 
ratio C/Mtot decreases monotonically with /i. Although it is as much as 30% 
greater than the An appearing in the statement of the hoop conjecture in the 
limit ^ 0, it is only about 2% less than An in the limit /i — > cx) (the limit of 
extreme prolateness). 

Does the absence of trapped surfaces on A/" when A < Xo{fi) suggest that 
trapped surfaces are completely absent from the spacetime? If so, this could be 
viewed as evidence that the singularity produced by the collapse of the shell is 
naked in this regime. However, we do not believe that the absence of trapped 
surfaces on Af when A < Ao(/i) provides strong evidence for the absence of 
trapped surfaces in the spacetime. When A = 1/8 and > 0, the outermost 
MTS is the intersection of Af with a hyperplane at a time t = ^mts- The 
intersection of A/ with any hyperplane with < < < tMTS is a trapped surface 
[0, so there is a large "trapped region" on Af. There is also a large trapped 
region to the future of Af. Namely, if T denotes a trapped surface on Af and £ 
denotes the null hypersurface generated by the outgoing null geodesies from T, 
then (by the Raychaudhuri equation) every cross-section of £ also is trapped. 

^In comparison, the "isoperimetric inequality" A < IGnM^^f. is satisfied by a wide margin 
for the MTS's occurring in the highly prolate case. Thus, the validity of this inequality [ho|| 
does not provide a "reason" why MTS's do not exist for A < Xo{fi). 
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As A is decreased (with fi held fixed) the trapped region on Af begins to "migrate 
away" from the singularity at t = near the poles {9,0' = 0). At A = Ai(/Lt), 
the trapped region on Af begins to migrate away from the singularity near the 
equator (z = 0) as well. For Xo{^) < A < Ai(/Lt), the trapped region on Af 
presumably is the region bounded by the two MTS's present on Af (see fig. ^). 
However, trapped surfaces also must continue to exist in the region to the future 
of swept out by outgoing null geodesies orthogonal to the trapped surfaces 
on Af. When A = Aq, there is a single MTS, A4, on A/", but no trapped surfaces 
on Af. However, since the shear of the outgoing null geodesies orthogonal to 
A4 is nonvanishing, there must continue to exist trapped surfaces to the future 
of Af. Thus, the disappearance of trapped surfaces on Af at X — Xq does 
not signal the disappearance of trapped surfaces from the spacetime, but rather 
their migration to the future oiAf. However, the issue of whether or not trapped 
surfaces continue to be present as A — > cannot be investigated until we have a 
much more detailed knowledge of the spacetime exterior to Af. 
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A Existence and uniqueness of solutions 

In this appendix, we outline a proof that in the BIL model, given any to > 0, 
there exists a unique value of m yielding an MTS on Af with /(O) = tg. (That 
is, we prove that the function m(t()) described in the text is well-defined for all 
^0 > 0.) The key ideas in the existence and uniqueness proof are due to David 
Ross. 

We fix M,L, and to, and leave m to be determined. We wish to find a 
solution of the following system of equations: 



f(0)^0 (38) 
/(O) = to (39) 
'sin0-(logy) =— (40) 



smOde V de^ ° ' } V 
, 2 ) - -^("^ 



VC^) - /(^) (41) 
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^„ogVO,|).-f,f, (42) 

(43) 

One can derive the following implicit solution for the first three equations 
(|7[|^) 

2qz 1 1 fw-l\ 

= + log (44) 



" = 1/7^ (45) 



where 



This solution generates certain values /(L/2) > and < df /dz{L/2) < 1; 
taking these values as given (determined by our choices of M, L, and Iq), we focus 
on (^-43). This system of equations, a 2nd-order ODE with three boundary 
conditions, would be overdetermined if we also specified m. However, we instead 
recast eqs. (^0[), (^2[), and ( ^3| ) in terms of the "scale-invariant" function V = 
V/m, which obeys 

^(») = (49) 

These equations do not involve m. The remaining boundary condition ( ^ ) is 
then satisfied if and only if m is chosen to be 

m = im. (50) 

It follows that if there exists a unique solution to (|47|-^9|), then there exists a 
unique value of m yielding an MTS with the chosen to. 

Now write k = —df /dz{L/2), let primes denote derivatives with respect to 
6, and define W = ([ogV)' . This gives 

^ {Wsiney^t (51) 



sin fc* V 

V' = WV (52) 
W{^) = k (53) 



16 



W(0) = (54) 

In order to prove the existence and uniqueness of solutions to this boundary- 
value problem, we convert it to an initial-value problem and use a "shooting 
argument", i.e. we consider solutions to ([5^), ( ]5^ ) and 

V{0) = Vo>0 (55) 

and then demonstrate that one can choose Vq uniquely in order to satisfy the 
boundary conditio n (p3|) . We will first prove local existence and uniqueness of 
solutions to (^2|), and (|5^) around 9 = 0, then show global existence 

and uniqueness on the interval [0, tt/2] for at least some choices of Vq, and finally 
argue that dsS) can always be satisfied. 



To begin, we note that IpWj and (54) are equivalent to the integral equation 




W = ^^l I 1 - -r^ ) sint (56) 



V{t) 

whereas eqs. (|52|) and (feSl) are equivalent to 



V{e) =Vo+ W{t)V{t) dt (57) 







To establish local existence of solutions to (^6|) and (q7|), we use a variant of 
"Picard iteration" . Consider the iteration scheme given by 



Wj{e) = ^- 1--^ \sintdt (58) 




V,{e) = Vo+ W,^i{t)Vj-i{t) dt (59) 





with initial "guesses" Vj=i = Vq and Wj^i = 0. It can be shown, by demon- 
strating that the Wj's and Vj's are bounded in some small interval beginning 
at 6* = 0, that the iterates converge uniformly as 

|V^.-VS-i|<^ (60) 

\W,-W,.,\<^^ (61) 

where C is a constant. Because the convergence is uniform, the limit functions 
V = limVj and W = limWj satisfy (^6|) and (^7|), and so constitute a (local) 
solution to the initial-value problem. 

Local uniqueness is proved similarly. Let Vi , W\ and V2 , ^2 both be solutions 
of the initial-value problem. Then, again using the fact that both of these 
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solutions are bounded on a sufficiently small interval starting at 6* = 0, one can 
show that for all positive integers j, 



\Vi-V2\<P^^^ (62) 

\m-W2\<P^^^p^ (63) 

where P and Q are constants. Therefore, Vi = Vz and Wi = W2, and solutions 
of the initial- value problem are unique near 9 = 0. 

By a standard theorem in the theory of ODE's then, this unique local 
solution to eqs. ( p6| ) and (|5^) can be extended away from 6 = until it becomes 
unbounded. It remains to be shown that Vq can be chosen so that the solution 
can be extended all the way to 6* = 7r/2 and satisfies VF(-|) = k. It is not 
difficult to prove that, for fixed 6 > 0,V and W are strictly increasing functions 
of Vq, so if we can find a Vq which yields the correct matching condition, this 
Vq will be unique. 

The analysis of global existence of solutions on the interval [0, 7r/2] naturally 
divides into the following three cases. 

• Vo > 4. Then V and W are strictly increasing functions of 9. Furthermore, 
V and W are bounded on the entire interval [0, 7r/2], according to 

Vo < V{9) < Vo exp ( / ^^7^ d?j (64) 

< W{9) < (65) 
smt^ 

with equalities occurring only at ^ = 0. As Vq 00, W{'k/2) 1, while 
as Vq ^ 4, W{-k/2) 0. By continuous dependence of solutions on initial 
data, it follows that all values of W{tt/2) satisfying < W{Tr/2) < 1 can 
be achieved. 

Vo = 4. Then V{9) = 4 and W{9) = 0. 

Vo < 4. Then < for 61 > 0, and, for fixed %, V{Vo,9) is a strictly 
decreasing function of 6*. (Recall that for fixed 9, ^(Vo, ^) is also a strictly 
increasing function of Vq.) In this case, it is possible that for some Vq, 
the solution "blows up" (F — > and W — > —00) before reaching 9 = tt/2, 
in which case a global solution does not exist. (It is not difficult to show 
that this is the only way in which a global solution can fail to exist.) 
We wish to show that given any k e (— oo,0), Vq can nevertheless be 
chosen to yield a global solution with W{tt/2) = k. First, we note that 
since = 4 is a global solution, it follows from continuous dependence of 
solutions on initial data that there exist global solutions on the interval 
[0, tt/2] for some Vo < 4. Let c be the greatest lower bound of the set of Vo 
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for which global solutions exist. Then, from the monotonicity properties 
of V, it follows that global solutions exist for all Vq € (c, 4). Let p be 
the greatest lower bound of V{n/2) for these global solutions (so that, 
by monotonicity, as Vq — > c, F(Vo7''i'/2) p)- Then there exist three 
possibilities: (i) li p — 0, then it can be shown that as F(Vb,7r/2) —^ 0, 
W{tt/2) — oo. Therefore, by continuous dependence of solutions on 
initial data, any value W{tt/2) £ (— oo, 0) can be obtained by some choice 
of Vo S (c, 4). (ii) If p > and c = 0, then we contradict the fact that 
V{9) is strictly decreasing in 9. (iii) li p > and c > 0, then one can 
show, using continuous dependence of solutions on initial conditions and 
the fact that solutions can fail to exist globally only by V approaching 
at some < tt/2, that there exists a global solution with Vq = c and 
V(tt/2) = p. However, again using continuous dependence of solutions on 
initial data, it then follows that there exist global solutions with Vq < c. 
Thus we obtain a contradiction. Consequently, only case (i) is possible. 

In summary, for any Vq > c, we can obtain a global solution to the initial- 
value problem on the interval [0, 7r/2], and we can uniquely choose Vq in this 
domain in order to achieve any W{tt/2) < 1. (Recall that we only needed 
— 1 < W{tt/2) < 0.) Therefore, there exists a unique solution to the boundary- 
value problem ( p7[p^ ), and consequently a unique value of m, determined by 
( ^0[ ) , yielding an MTS with the chosen ■ 



B Limiting behavior of m{to) 

Given M and L, we showed in appendix A that m is determined by a choice of 
to- This dependence is mediated by the "matching conditions" between / and 
V, the center and end solutions. We write 

Xito)= f{to,z = L/2) (66) 

df 

Y{to)^-^{to,z = L/2) (67) 

where /(<o, z) denotes the center solution, f{z) (see eqs. (^),(^), with /(O) = 
to- It is convenient to view X and Y as independent variables and to view m 
as a function of X and Y, 

m^m{X,Y) (68) 

We may then write 

dm dm dX dm dY . , 

dt^^dXdt^^dYdt^ ^ ' 

In this appendix, we derive the limiting behavior as to ^ a-nd as tg — > oo 
of each of the terms on the RHS of this equation, and so discover the limiting 
behavior of m(to). 

We first examine how m changes with X when Y is held fixed. Since m is 



determined by solving the system of equations (^-49), and these equations are 
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independent of X, it follows that if X is changed, keeping Y fixed, we obtain 
the same solution V. From cq. (dOI), we therefore obtain 



dm 1 

Next, we change Y while holding X fixed. Since we found in the previous 
appendix that V{Tr/2) and W{Tr/2) are strictly monotonically increasing func- 
tions of Vq, an increase in Y results in a decrease in W{7t/2), a decrease in Vq, 
and a decrease in V{n/2). By eq. (^o|), the required m must increase, and so 

(ri) 

For the remaining terms, we need to refer to the implicit solution (Q) and 
(Esl) for f{z), as well as the first integral of eq. (pTh 



By taking the derivatives with respect to tg of the implicit solution (Q) and the 
first integral d72) we obtain, respectively. 



dto to 2ta\l X ^ ' 

= (74) 

Although our knowledge of dm /dto is limited by the fact that we have only 
qualitative knowledge of dm/dY, eq. (|7|), we can nevertheless derive its exact 
behavior in the limits tg and Iq oo. When tQ <C i, the implicit solution 
for / yields, to lowest order, 

X.f.|,o.(i|) (75, 

implying 

dtf) \ aL 

On the other hand, since < df /dz{to, z) < a, 

to< X <to + aL/2 (77) 

Then eq. (|7|) gives 

dto 



log ("iL ) ^ _oo as io ^ (76) 



dX 

1 as to oo (78) 
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Since, by the above, X aL/2 as io ^ and X ^ oo as ^ ^ oo,we have 
from eq. (Q) 

dY / -1 as to ^ ^^g^ 



dtQ \ as ^0 ^ oo 

According to eq. ([72|), 



[ as to ^ oo ^ ' 

Now, applying a result of the previous appendix, y(7r/2) approaches some value 
between and 4 in the first limit, and y(7r/2) goes to 4 in the second limit. 
Thus 

dm ( const. > i as to ^ 
dX \ J as to oo 

Returning to eq. (|69|), we have that 

> -oo as to ^ (82) 

"to 

whereas m approaches a finite limiting value, mi, as to — > (since there are no 



(81) 



special difficulties encountered in solving (47-49|) when df /dz{L/2) = a). In 
the other limit, we have simply 

dm 1 

_^_asto-oo (83) 

Finally, note that when m = 0, we can explicitly solve for V{6), finding that 
(log F)'(7r/2) = 1. However, it is impossible to choose to in order to satisfy 
the corresponding matching condition df /dz{L/2) ~ —1. Thus there exists no 
MTS confined to the null hypersurface Af when m = 0. 

The above results imply that m(to) has the behavior described in section 4 
(see figure || above). 
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